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Several statistical mechanical theories predict that colloidal suspensions of highly charged 
macroions and monovalent microions can exhibit unusual thermodynamic phase behavior when 
strongly deionized. Density-functional, extended Debye-Hiickel, and response theories, within mean- 
field and linearization approximations, predict a spinodal phase instability of charged colloids below 
a critical salt concentration. Poisson-Boltzmann cell model studies of suspensions in Donnan equi- 
librium with a salt reservoir demonstrate that effective interactions and osmotic pressures predicted 
by such theories can be sensitive to the choice of reference system, e.g., whether the microion density 
profiles are expanded about the average potential of the suspension or about the reservoir potential. 
By unifying Poisson-Boltzmann and response theories within a common perturbative framework, it 
is shown here that the choice of reference system is dictated by the constraint of global electroneu- 
trality. On this basis, bulk suspensions are best modeled by density-dependent effective interactions 
derived from a closed reference system in which the counterions are confined to the same volume as 
the macroions. Lower-dimensional systems (e.g., monolayers, clusters), depending on the strength of 
macroion-counterion correlations, may be governed instead by density-independent effective interac- 
tions tied to an open reference system with counterions dispersed throughout the reservoir, possibly 
explaining observed structural crossover in colloidal monolayers and anomalous metastability of 
colloidal crystallites. 



I. INTRODUCTION 

A variety of experiments have demonstrated 
the unusual thermodynamic properties of deionized 
charge-stabilized colloidal suspensions [1]. Aqueous 
suspensions of highly charged macroions and mono- 
valent microions at sub-millimolar ionic strengths 
reportedly can display liquid- vapor coexistence [2], 
stable void structures [3-5], compressed crystals [5, 
6], metastable crystallites [7], and macroion gath- 
ering near glass plates [7, 8]. Such phenomena 
have been interpreted by some workers [2-5] as evi- 
dence for pair attraction of like-charged macroions, 
in seeming defiance of the classic Derjaguin- 
Landau-Verwey-Overbeck (DLVO) [9] and Poisson- 
Boltzmann (PB) [10] theories. Others have at- 
tributed anomalous behavior to nonequilibrium phe- 
nomena [11], polyelectrolyte impurities [12], and 
many-body effects [13]. While some molecular sim- 
ulations of the primitive model display macroion ag- 
gregation [14-23], such computationally demanding 
methods are usually limited to size and charge asym- 
metries corresponding to relatively strongly corre- 
lated microions and weakly charged macroions. 

Several common statistical mechanical theo- 
ries, including density- functional [24-27], extended 
Debye-Hiickel [28-30], and response [31-35] theories, 
predict similarly surprising thermodynamic phase 
behavior of deionized suspensions [36] . Such coarse- 
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grained theories preaverage the microion degrees of 
freedom, reducing the ion mixture to an effective 
one-component model. Most practical implemen- 
tations assume a mean-field approximation for the 
microion structure, which neglects microion corre- 
lations; some form of linearization approximation, 
which ignores nonlinear screening of macroions by 
microions; and a fixed (state-independent) effective 
macroion charge. Under these assumptions, the var- 
ious approaches reduce to variants of linearized PB 
theory, all predicting a screened-Coulomb (Yukawa) 
effective pair potential and a one-body volume en- 
ergy. Although independent of macroion coor- 
dinates, the volume energy contributes density- 
dependent terms to the total free energy that can 
drive a spinodal instability of highly charged suspen- 
sions below a critical salt concentration [24-28, 35]. 
Such unusual phase behavior occurs, however, at pa- 
rameters as yet inaccessible to primitive model sim- 
ulations and may be qualitatively modified by ion 
correlations and nonlinear screening. 

Recent studies of colloidal suspensions in Don- 
nan equilibrium [37, 38] with an electrolyte reser- 
voir demonstrate that predictions of linearized PB 
theory can be sensitive to the choice of reference 
system [39-42]. While expansion of the microion 
densities about the average potential of the system 
leads to phase separation at low salt concentrations, 
expansion about the reservoir potential strictly pre- 
dicts phase stability, in qualitative agreement with 
nonlinear PB theory [39-43]. These studies suggest 
that the predicted phase separation may be merely 
a spurious artifact of linearization approximations. 
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Although nonlinear PB theory is often presumed 
superior to its linearized form and, by implication, 
to related linearized theories, it is well to remem- 
ber that mean-field theories do not faithfully model 
nonlinear screening near highly charged macroions, 
where counterions are strongly correlated. Indeed, 
strong countcrion association can rcnormalizc the 
effective macroion charge [44-50], modifying inter- 
particle interactions and thermodynamic properties. 
Furthermore, nonlinear PB theory is computation- 
ally practical only within cell models [51], which re- 
quire computing the microion distribution around 
just a single macroion, but thereby completely ne- 
glect macroion correlations. In contrast, linearized 
theories predict analytical expressions for effective 
interactions that are independent of any artificial cell 
geometry, become accurate sufficiently far from the 
macroions, and can be input into separate theories 
or simulations that incorporate ion correlations and 
charge renormalization. Linearized theories thus 
may play a vital role in multiscale approaches to 
modeling charged colloids. 

The main purpose of this work is to carefully 
analyze implementations of linearization approxi- 
mations in coarse-grained theories of charged col- 
loids. Connections between PB and response the- 
ories are established and exploited to show that 
the constraint of global electroncutrality dictates 
the optimal choice of reference system. Theories 
linearized about a closed reference system, whose 
counterions are confined to the same volume as the 
macroions, are shown to predict density-dependent 
effective interactions and - assuming a fixed effec- 
tive macroion charge - phase instability of deion- 
ized suspensions. In contrast, linearization about 
an open reference system, with counterions dis- 
persed throughout the reservoir, predicts density- 
independent effective interactions and phase stabil- 
ity, but violates electroneutrality for bulk suspen- 
sions. Lower-dimensional suspensions (e.g., mono- 
layers, clusters), on the other hand, may be modeled 
best - for sufficiently weak macroion-counterion cor- 
relations - by an open reference system, with impor- 
tant implications for thermodynamic properties. 

The remainder of the paper is organized as fol- 
lows. Section II first defines the primitive and ef- 
fective one-component models. In Sec. Ill, effective 
electrostatic interactions and osmotic pressures are 
then derived from response theory and PB theory 
for closed and open reference systems. Correspond- 
ing predictions for the equation of state and phase 
diagram are presented and contrasted in Sec. IV. 
Section V closes with a summary and conclusions. 




FIG. 1: Colloidal suspension (left) of macroions (larger 
circles) and microions (smaller circles) in Donnan equi- 
librium with a microion reservoir (right) via a semiper- 
meable membrane (dashed line), which allows exchange 
only of microions. 



II. MODELS 



Primitive Model 



We consider a charge-stabilized colloidal suspen- 
sion in Donnan equilibrium [37, 38] with a microion 
reservoir (e.g., electrolyte). As depicted in Fig. 1, 
the suspension and reservoir are separated by a 
semipermeable membrane, which allows exchange of 
solvent and microions, but not macroions. The vol- 
umes of the suspension and reservoir are denoted 
by V and V r , respectively. Within the primitive 
model of charged colloids [1], the solvent is mod- 
eled as a dielectric continuum, of dielectric constant 
e, the macroions as N m charged hard spheres, of ra- 
dius a and effective valence Z (charge — Ze) and the 
microions as monovalent point charges, N + positive 
and N_ negative, totaling JVp = N + + N^ in the sus- 
pension. It is assumed that all ions interact via only 
bare Coulomb and cxcluded-volume pair potentials 
and have the same dielectric constant as the solvent, 
justifying neglect of polarization effects. 

The macroions have an average number density 
n m = N m /V, their hard cores occupying a volume 
fraction rj = (47r/3)n m a 3 . The microions have av- 
erage number densities n± = N±/V in the free 
volume V = V(l — rj) outside the macroion cores. 
The reservoir is assumed to be a 1:1 symmetric 
electrolyte with number density of salt ion pairs 
n r . Given Z counterions per macroion and N s salt 
ion pairs, the ion numbers are related according to 
N + = ZN m + N s and AT_ = N s and the ion densities 
according to n+ = Zn m /(1 — ?/) + n s and n_ = n s , 
where n s = N B /V is the average salt density in the 
free volume. Overall electroneutrality of the suspen- 
sion imposes the constraints N + — A_ = ZN m and 
n + - n_ = Zn m /(1 - rj). 
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B. Effective One-Component Model 

By averaging over the microion degrees of free- 
dom, the primitive model is mapped onto a one- 
component model of pscudomacroions, governed by 
effective interactions [52-56]. This coarse-graining 
procedure acts on the semigrand partition function 
((exp(— f3H))fj,) m , where ( ) m denotes a canonical 
trace over macroion (m) coordinates, ( ) a grand 
canonical trace over microion (/i) coordinates, H 
is the Hamiltonian, and (3 = 1/ksT at tempera- 
ture T. Splitting H into macroion, microion, and 
macroion-microion interaction terms, according to 
H = H m + + H m+ + H m _, and tracing over the 
microions yields the semigrand potential, 

{l aB = -k B Tln(exp(-0H eB )) m , (1) 

where H c s = H m + fi M is an effective one-component 
Hamiltonian and 

= -k B Tln (cxp [-PiHp + H m+ + 

(2) 

is the grand potential of the microions amidst fixed 
macroions. Equations (1) and (2) provide a for- 
mal and exact foundation for both response theory 
(Sec. Ill A), which is based on a perturbative expan- 
sion of flfi about a reference system, and PB theory 
(Sec. IIIB), which is based on a density- functional 
expansion of the grand potential functional. 

III. THEORY 

A. Linear-Response Theory 

Response theory of charged colloids [31-35] views 
the electrostatic potential of the macroions as an 
"external" potential that perturbs the microions, in- 
ducing structure in their (otherwise uniform) den- 
sity profiles. Within this view, the microion grand 
potential [Eq. (1)] can be expressed as 

n^ = n + f d\ ((H m+ ) x + (# m _> A ) , (3) 

Jo 

where flo = — fcsTln (exp(— PH^)) is the grand po- 
tential of a reference system in which the macroions 
are uncharged and the microions are unperturbed, 
A is a charging parameter, and ( ) x represents a 
grand canonical ensemble average for macroion va- 
lence XZ. In practice, it proves convenient to add to 
and subtract from Qq the energy Ej, of a uniform, 
neutralizing background, having charge density op- 
posite that of the unperturbed microions, and to de- 
fine flp = ilo+Ei, as the grand potential of a uniform 
(electroneutral) reference microion plasma. 




FIG. 2: Open (left) and closed (right) reference systems 
in linear-response theory of a salt-free colloidal suspen- 
sion. Larger circles represent uncharged (hard-sphere) 
colloids and smaller circles charged counterions. 

For a sealed suspension whose macroions and mi- 
croions share the same volume, the reference system 
is unambiguous. For a suspension in Donnan equi- 
librium, however, two choices seem to be possible: 

(1) A "closed" reference system, in which the coun- 
terions are confined, with the macroions, to the sus- 
pension. Independent of reservoir size, the microion 
densities in the reference system are then equal to 
those in the suspension. 

(2) An "open" reference system, in which the mi- 
croions arc uniformly distributed throughout the 
combined volume of the system and the reservoir. 
In the limit of an infinite reservoir (V r /V — > oo), 
the microion densities in the reference system equal 
those in the reservoir. 

The two choices are depicted in Fig. 2 in the ide- 
alized case of zero salt concentration. The second 
choice may seem most natural, since in the presence 
of neutral macroions the counterions would diffuse 
throughout the available volume to maximize en- 
tropy. On the other hand, in the real system (with 
charged macroions), electroncutrality confines the 
counterions to the suspension. Does it matter which 
reference system is adopted? By considering each 
case in turn, it is shown below that the choice mat- 
ters, in practice, only for strongly deionized suspen- 
sions, but that the optimal choice then may depend 
on system dimensionality. 

Linear-response theory [31-35] is first reviewed, 
emphasizing the role of the reference system. The 
external potential at position r is defined as 

Uext(r) = / dr'w m± (|r-r'|)n m (r'), (4) 
Jv 

where v rn ±(r) = Ze 2 /er (r > a) are the macroion- 
microion Coulomb pair potentials and n m (r) is the 
macroion density. The macroion-microion interac- 
tions then can be expressed as 

H m ± = / drw oxt (r)n ± (r) 
Jv 

= y;^v m ±(k)fi m (k)n±(-'k), (5) 

k 
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where n±(r) denote the microion density profiles 
and the Fourier transform is defined according to 
n TO (k) = J v , dr cxp(— ik ■ r)n m (r), with inverse 

n m (r) = (l/V)Z)k ex P( ik ' r)"m(k). 

Expanding n±(r) in a functional Taylor series in 
powers of v cx t (r) yields [33] 

h±(k) = X ±(k)v m ±(k)n m Qc) + ■ ■ ■ , fc^O, (6) 

where x±(fc) are (partial) linear-response functions, 
associated with the microion structure in the ref- 
erence system, higher-order terms involve nonlin- 
ear response functions and higher powers of the 
macroion density [34, 35], and the zero-wavevector 
(long- wavelength) component is fixed by the normal- 
ization condition n±(0) = N±. Combining Eqs. (3), 
(5), and (6), the effective Hamiltonian can be writ- 
ten as [32-34] 



H cf f = Hn S + E + - ^2 v eS (rij) 



(7) 



where -ffHS is the hard-sphere (HS) Hamiltonian as- 
sociated with the macroion cores, E is a one-body 
volume energy, v e s(r) is an effective electrostatic 
pair potential between macroions whose centers are 
separated by r, and higher-order terms involve effec- 
tive many-body interactions [34, 57, 58]. The effec- 
tive pair potential, 

v e s{r) = v mm (r) + v ind (r) (8) 

is the sum of the bare macroion-macroion Coulomb 
pair potential v mm (r) = Z 2 e 2 /er and a microion- 
induced pair potential v ind (r). Within a linear- 
response approximation [33], 

i>ind(k) = x (k)[v m+ (k)} 2 , (9) 

where x(k) = X+(k) + X-(k) is the (total) linear- 
response function. Equations (7)-(9) are general ex- 
pressions for the effective interactions. More explicit 
expressions depend on the specific reference system. 

In a closed reference system, where the unper- 
turbed microions have average densities n± , a simple 
ideal-gas approximation yields 













m 


- 1 




ln ( n r ) 


- f 













(10) 

assuming weakly correlated microions with (elec- 
trochemical potentials [59] fixed by the reservoir: 
M± = /'r = fcsTln(n r A 3 ), A being the microion 
thermal wavelength. The one-body volume energy 
is given by [33] 

E = a p + ^ind(O) + N m (n+ - »_) 



x lim 

fe^O 



1 Z 
v m +{k) - — Wi„d(fc) + -v{k) 



(11) 



where v(k) is the Fourier transform of the 
counterion-counterion Coulomb pair potential 
v(r) = e 2 /er. The linear-response functions are 
related to the static structure factor S(k) of the 
closed reference microion plasma [60] via 



X±(k) = TPn±S(k) 



1 - n^c(k) ' 



(12) 



where = n + + n_ = Zn m /(1 — rf) + 2n s is 
the total average microion density in the free vol- 
ume of the suspension and c(k) is the Fourier trans- 
form of the counterion-counterion direct correla- 
tion function. Within a random-phase approxima- 
tion [60] for the microion structure, which neglects 
all but asymptotically long-ranged microion correla- 
tions, c(k) = —f3v{k) and the linear-response func- 
tions become 



X±(k) 



(3n± 



l + K 2 /k 2 



(13) 



where k = ^Attu^Xb is the Debye screening con- 
stant and \b = /3e 2 /e is the Bjerrum length. Tak- 
ing into account exclusion of the microions from the 
macroion cores leads to macroion-microion poten- 
tials of the form [33] 



(k) = T 



4nZe 2 



e(l + na)k 2 



cos(fca) 



sin(fca) 



(14) 

Substituting Eqs. (13) and (14) into Eqs. (9) and 
(11) results in explicit analytical expressions for the 
effective interactions, namely a screened-Coulomb 
(Yukawa) effective pair potential [33] 



v e f[(r) 



Z 2 e 2 { e Ka \ 2 e- Kr 



e \ 1 + na 



r>2a (15) 



and a volume energy 



(3E = (3fl f 



N m Z 2 k\ 



B 



N m Z n A 



1 + na 



; (16) 



whose three terms incorporate, respectively, the mi- 
croion entropy, the macroion self-energy, and the av- 
erage microion potential energy. The latter term is 
a direct manifestation of the Donnan effect [37, 38] 
- the unequal distribution of microions between the 
suspension and reservoir resulting from the imper- 
meability of the interface to the macroions. Dif- 
fusion of counterions from the suspension into the 
reservoir generates an interfacial charge and a re- 
sultant electric field that pulls the counterions back 
to maintain global electroneutrality of the suspen- 
sion. The electrostatic potential of the suspension is 
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thereby shifted, by the Donnan potential typ, rela- 
tive to that of the reservoir to equalize the microion 
chemical potentials. In terms of ^n, the final term 
in Eq. (16) can be expressed as (n + — n_)e x l , D/2, 
which is simply the work required to move microions 
from the reservoir (at zero potential) to the suspen- 
sion (at potential ^>d)- 

It is vital to note that the effective pair potential 
[Eq. (15)] and volume energy [Eq. (16)] derived from 
the closed reference system depend implicitly (via k) 
on the average microion densities in the suspension. 
The electroncutrality constraint then imposes a de- 
pendence of Vett(r) and E on the macroion density. 
As discussed below in Sec. IV, this density depen- 
dence entails an effective many-body cohesion that 
can profoundly impact the bulk phase behavior of 
countcrion-dominatcd suspensions. 

In an open reference system, the unperturbed mi- 
croions are not confined with the macroions, but 
instead uniformly fill the combined free volume of 
the suspension and reservoir. The effective inter- 
actions are now still described by Eqs. (15) and 
(16), but with the average microion densities in 
the suspension n± replaced by the average mi- 
croion densities in the suspension + reservoir, i.e., 
(n± + n r V r /V')/{l + V r /V). In the limit of an infi- 
nite reservoir (V r /V' — > oo), the reference microion 
densities approach n r and the reference microion 
plasma grand potential [Eq. (10)] j3Vl p — > —2n r V. 
In this same limit, the linear-response functions, now 
associated with the microion structure of the reser- 
voir, become [cf. Eq. (13)] 



{3n r 



1 



l/k 2 



(17) 



where n r = ^J&nn r Xb is the reservoir screening con- 
stant. Likewise, the effective pair potential retains 
the Yukawa form [Eq. (15)] with n replaced by n r 
and the volume energy reduces to [cf. Eq. (16)] 



(3E r = -2n r V' 



D 



1 



(18) 



Note that the electroneutrality term [final term of 
Eq. (16)] vanishes, since in the infinite-reservoir limit 
the reference microion plasma contains equal densi- 
ties of positive and negative microions, the countc- 
rions having "evaporated" into the reservoir. 

Thus, when applied to a suspension in Don- 
nan equilibrium with a microion reservoir, linear- 
response theory predicts nontrivial dependence of 
the effective interactions on macroion density only 
when the effective Hamiltonian is perturbed about a 
closed reference system. If the perturbative expan- 
sion is performed about an open reference system, 



the volume energy - in the case of an infinite reser- 
voir - is independent of macroion density and hence 
does not influence thermodynamics. Clearly the dis- 
tinction between reference systems can have practi- 
cal relevance only if the counterion concentration is 
not overwhelmed by the salt ion concentration. Sec- 
tion IV identifies a physical criterion for selecting the 
optimal reference system and explores implications 
for phase stability of deionized suspensions. 



B. Linearized Poisson-Boltzmann Theory 

An alternative, and widely studied, approach to 
modeling charged colloidal suspensions and polyelcc- 
trolyte solutions is the Poisson-Boltzmann theory, 
here briefly reviewed with emphasis on connections 
to response theory and the role of the reference sys- 
tem. Poisson-Boltzmann theory can be derived from 
the microion grand potential [61] 

n M [n+(r),n_(r)] =T» - fi+N+ - M _V_, (19) 

regarded as a functional of the nonuniform microion 
number densities n±(r) in the external potential 
Vext(r) of the macroions [Eq. (4)]. The microion 
Helmholtz free energy functional separates natu- 
rally, according to 

.7>[n + (r),ra_(r)] = T\ A + T cxt + F ex , (20) 

into an ideal-gas free energy 

■Fid = k B T [ dr (n+(r){ln[n+(r)A 3 ] - 1} 
Jv 

+ n_(r){ln[n_(r)A 3 ]-l}), (21) 
due to the microion entropy, an external free energy 



dru oxt (r)[n + (r) - rc_(r)], 



(22) 



associated with microion-macroion interactions, and 
an excess free energy T ex , due to microion- microion 
interactions. In a mean-field (Hartree) approxima- 
tion, which neglects all interparticlc correlations, the 
excess free energy can be written as 

F C x = - f dr / dr'vflr - r'|)[n+(r) -n_(r)] 
* Jv Jv 

x [n + (r')-n_(r')]. (23) 

Combining Eqs. (19)-(23) and minimizing 
f2 M [n+(r),n_(r)] with respect to n±(r), for a given 
macroion density, leads to Boltzmann distributions 
for the equilibrium microion densities: 



n±(r) = n r exp[=F/3e^(r)], 



(24) 



G 



where the electrostatic potential ^(r) is defined via 



e*(r) 



dr'u(|r-r'|)n(r'), 



(25) 



with n(r) = n + (r) — ?i_(r) — «/(r) the total num- 
ber density of all charges, including those fixed on 
the macroion surfaces n/ (r) . Equation (24) and the 
exact Poisson equation 



V 2 *(r) 



Aire 



n(r) 



(26) 



together determine 4 r (r) and n±(r) under prescribed 
boundary conditions. The total free energy func- 
tional (for fixed macroions), T = E mm + JF^, which 
includes the macroion-macroion Coulomb interac- 
tion energy, E mm = J2kj v mm(nj), is given by 

T = Tid 



dr ^(r)n(r) 



or, within a cell model approximation, 



T = T id + N rn — / dr|W(r 

™ ./cell 



(27) 



(28) 



using Eq. (26) and assuming vanishing electric field 
(V\& = 0) at the cell boundary. 

More generally, independent of any cell geome- 
try, and for any fixed macroion distribution, the mi- 
croion free energy functional can be expressed as 

= F id + ^^v m+ (k)h m (k)[h + (-k) - h_(-k)} 

k 

+ ^E {i ( fc )["+(k) - Mk)][n+(-k) - n-(-k)]. 

(29) 

A systematic scheme for further approximations is 
based on expansion of the microion density profiles 
[Eq. (24)] about a reference potential V&o and a func- 
tional Taylor-series expansion of the ideal-gas free 
energy [Eq. (21)] about reference densities n±\ In 
linearized PB theory, the expansions of n±(r) are 
truncated at linear order in the potential deviation, 



n ± (r)=4 0) {lT/3e[*(r)-* ]}, 



(30) 



and the expansion of .Fid is truncated at quadratic 
order in the density deviations, 

F id = ksT^ jiV t ln(n| 0) A 3 )-^ 0) 
i=± I 



2n (0) JV 



dr 



n, (r) - n 



(0) 



(31) 



where = n r exp(=p/3e v E'o) ar< 5 the reference mi- 
croion densities. Analogous to the choice in linear- 
response theory between closed and open reference 
systems, here two choices of reference potential and 
corresponding microion densities are apparent: 

(1) The average potential and microion densities in 
the suspension, "fo = W and n± = n±. 

(2) The average potential and microion densities in 

the reservoir, $o — and n± = n r . 

Although the second choice is most commonly as- 
sumed, the first has been advocated [39-41] as 
better exploiting the supposedly weak deviations, 
[^(r) - ¥] and [n±(r) - n±]. 

Expanding n± (r) first about "to = an d combin- 
ing the Fourier transform of the linearized expansion 
with that of Eq. (25), 

e$(k) = v m+ (k)n m (k)+v(k)[n + (k)-n-(k)], (32) 

the microion densities take on the form 

n±(k) = =Fn±/3e#(k) = x ±(k)v m ±(k)h m (k), k ^ 0, 

(33) 

where x±(&) arc identical to the linear-response 
functions that appear in Eqs. (6) and (13). Equa- 
tion (33) implies that 

n + (k) - n_(k) = x(k)v m +(k)h m (k), k ^ 0. 

(34) 

From the inverse transform of Eq. (33), using 
Eqs. (13) and (14), the electrostatic potential around 
an isolated macroion is then given by 



tf(r) = -■ 



t 1 + Kd r 



r > a, 



(35) 



which is identical to the solution of the Pois- 
son equation [Eq. (26)] with boundary conditions 
#'(r)lr=a = -Ze/ea 2 and (r) -> as r -> oo. The 
corresponding expansion of Fid about the average 
microion densities becomes 



Fi d = k EsT^2 \ Ni[ln(mA 3 ) - 1] 

i=± 



k#0 J 



(36) 



where the k = term vanishes via normalization: 
f v ,drn±(r) = n±V. Combining Eqs. (29), (33), 
and (36) then yields the microion Hchnholtz free en- 
ergy (equilibrium value of .F M ) for fixed macroions, 
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7^77 ^2 v ind (k)h m (k)h m (-k) + (n+ - n_) 



lira 



N-)v(k) +N m v m +(k) 



(37) 



the fc = terms being determined by the elec- 
troneutrality constraint. Significantly, this free en- 
ergy generates precisely the same effective Hamilto- 
nian, H c h = _ffns + F — N^[i r , as derived above 
(Sec. Ill A) from response theory perturbed about a 
closed reference system [Eqs. (7)-(ll)]. 

Expanding n±(r) now about the reservoir poten- 
tial (^o = 0), the Fourier transform of the linearized 
microion densities combined with Eq. (32) yields 



n±(k) = =pn r ^e*(k) 



X ( ±\k)v^ ± (k)h m (k), fc^O, 
(38) 

from which 

n+(k) - fi_(k) = X {r) (k)v^ + (k)h m (k), k^0 7 

(39) 

where x± (&) are the linear-response functions of 

Eq. (17), X M (fc) = xi r) (fc) - X M (fc), and C^(fc) 
are given by Eq. (14), with k replaced by n r . The 
expansion of .Fid about the reservoir microion den- 
sity becomes 

A^[ln(ra r A 3 ) - 1] -n r V 



2n r V 



^J^fli^-k), (40) 

=± k 

the k = term now nonvanishing. Substituting 
Eqs. (38)-(40) into Eq. (29) then yields 

A^[ln(n r A 3 ) - 1] 



0Fn 



(N + + N_y 

UrV + 4n r V 



2yJ X]^nd( fe ) ft m(k)n m (-k), 
k 



(41) 



where 



( r ^(fc)[w^(fc)] 2 is an induced po- 



tential that depends (via K r ) on only the reser- 
voir microion density n r . This free energy corre- 
sponds to essentially the same effective Hamiltonian 
as derived from response theory perturbed about 
the open (reservoir) reference system, including the 
same effective pair potential [Eq. (15) with reservoir 
screening constant n r ], but a slightly different vol- 
ume energy [cf. Eq. (18)], 



0E r = V 



n,, - rir - 



N m Z 2 n r \i 



1 + K r a 



Two main conclusions follow from the above 
derivations. First, linearization about the aver- 
age potential and microion densities of the sus- 
pension leads to effective interactions that de- 
pend on macroion density, while linearization about 
the reservoir yields density-independent interac- 
tions. Second, linear-response and linearized PB 
theories are formally equivalent [62], with direct 
correspondences between closed/open and suspen- 
sion/reservoir reference systems. In Sec. IV, these 
connections are exploited to bolster previous argu- 
ments [39-42] for linearizing PB theory about the 
average potential of the suspension, rather than that 
of the reservoir. 



C. Thermodynamic Properties 

Despite the presence of three ion species in the real 
suspension, thermodynamic properties of the effec- 
tive one-component model depend on the chemical 
potentials of only pseudomacroions [L m and salt ion 
pairs ix s , since electroneutrality permits exchange of 
ions only in electroneutral units. A suspension in 
Donnan equilibrium with a salt reservoir, at fixed 
salt chemical potential fi s — 2/i,. = 2/csTln(n r A 3 ), 
has a total pressure p = n m fi m — to and a pseu- 
domacroion chemical potential \i m = (du> / dn m ) Ug , 
where lu = Ct/V. The semigrand potential density 
of the suspension can be expressed as 



w(n m ,(j, s ) = f e s + ■ 



(43) 



where f c g is the free energy density of the one- 
component system interacting via the Yukawa ef- 
fective pair potential [Eq. (15)], with screening con- 
stant k or k t , and e = E/V [Eq. (16)] or E r /V 
[Eq. (18)] for closed or open reference systems, re- 
spectively. The salt density in the suspension is 
determined by solving for n s the implicit relation 
fi s = (df/dn s )n m , where 



f{n m ,n s ) =u3 + n a n a 



(44) 



is the total free energy density of the suspension [63] . 

The dependence of effective interactions on the 
choice of reference system has important implica- 
tions for the osmotic pressure of deionized suspen- 
sions. Substituting the appropriate volume energies 
into Eq. (43) yields, for the closed reference system, 

n r, i \ Z(n + — n^)n\ B 

pp = n fl + n m + fip ex .{n ll ) — — ^ — , (45) 

4(1 + nay 

and for the open reference system and an infinite 
reservoir, 



(42) 



Pp = 2n r + n m + /3p cx (n r ), 



(46) 
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where 



IV. RESULTS AND DISCUSSION 



fdfe 



(47) 



N B /N n 



is the excess pressure due to effective macroion- 
macroion pair interactions (with respective screen- 
ing constant k or n r ) and / cx = F CK /V is the ex- 
cess free energy density. The first two terms on the 
right sides of Eqs. (45) and (46) are the pressure 
contributions from microion and macroion transla- 
tional entropy. The final term in Eq. (45), no- 
tably absent from Eq. (46), results from the den- 
sity dependence of the microion-macroion interac- 
tions and thus is a manifestation of electroneutral- 
ity. For the open reference system and a finite reser- 
voir, Eq. (45) applies with the simple reassignment 
n± — > (n± + n r V r /V')/(l + V r /V), i.e., average mi- 
croion densities in suspension + reservoir. 

The excess free energy density can be accurately 
approximated by a variational method [24-26, 35] 
based on first-order thermodynamic perturbation 
theory with a hard-sphere reference system [60]: 



fcx{n m ,n s ) = min <! / H s (n m , n s ; d) + 2m%\ 
(d) 



x J drr gus(r,n m ;d)v eS (r,n m ,n s )> , (48) 

where the effective hard-sphere diameter d is 
the variational parameter and /hs ( n m ,n s ;d) and 
<?Hs( r ) n m] d) are the excess free energy density and 
(radial) pair distribution function, respectively, of 
the HS fluid, computed here from the essentially 
exact Carnahan-Starling and Verlet-Weis expres- 
sions [60]. From the Gibbs-Bogoliubov inequal- 
ity [60], minimization with respect to d generates 
a least upper bound to the free energy. 

The phase diagram can be computed from a 
(Maxwell) common-tangent construction on the 
curve of lo vs n m at fixed salt chemical potential, 
which imposes equality of the pressure and of the 
chemical potentials of macroions and salt in coex- 
isting phases. As an internal consistency check, 
it is possible also to calculate the chemical po- 
tentials of the individual microion species, /x± = 
(dF/dN±)y,N m ,N T - Chemical equilibrium between 
the suspension and reservoir requires 



^ + - M _)=ln(^±)-2^+-i = 0, (49) 
n_ / n,. 



which follows from Eq. (16) and symmetry of k with 
respect to interchange of n + and n_ . 



The importance of choosing an appropriate ref- 
erence system in linearized theories becomes clear 
upon examining the osmotic pressure II = p — 
2n r ltBT - the difference between the suspension 
and reservoir pressures - at constant salt activ- 
ity z s = exp(/3/x s )/A 3 , along with the correspond- 
ing phase behavior. Figure 3 shows sample predic- 
tions of linear-response theory, with both closed and 
open reference systems, for macroion radius a = 50 
nm and fixed effective valence Z = 500. While 
the present theory neglects any influence of charge 
renormalization on Z, and thus on phase behavior, 
the results at least qualitatively illustrate the signif- 
icance of the clcctroncutrality constraint for deion- 
ized suspensions. 

Linearization about a closed reference system pre- 
dicts, at sufficiently low salt activity, a van der Waals 
loop in the equation of state (II vs if), i.e., a range 
of volume fractions over which the compressibility is 
negative, independent of reservoir size [solid curve 
of Fig. 3(a)]. This unusual spinodal instability, first 
discovered by van Roij et al. [24-26], implies sep- 
aration into macroion-rich (liquid) and macroion- 
poor (vapor) bulk phases [67] within the binodal 
of the corresponding phase diagram [solid curve of 
Fig. 3(b)]. Close examination reveals the instability 
to be a many-body effect driven by the density de- 
pendence of both the macroion self-energy and the 
electroneutrality terms in the volume energy [last 
two terms of Eq. (16)]. 

Linearization about an open reference system pre- 
dicts, in sharp contrast, a van der Waals loop that 
narrows, and a liquid-vapor binodal that shrinks, 
with increasing reservoir volume [dashed curves of 
Fig. 3]. For V r > V, the loop closes entirely and 
the binodal collapses, the compressibility then be- 
ing strictly positive at all salt concentrations, imply- 
ing a single stable fluid phase. Such a trend would 
imply the unphysical possibility of influencing bulk 
phase behavior by controlling the exchange of mi- 
croions between suspension and reservoir, i.e., sim- 
ply adjusting chemical boundary conditions. The 
corresponding salt concentrations in the suspension 
c s (inset to Fig. 3(a)) also show qualitatively distinct 
trends, decreasing or increasing with increasing r\ for 
closed or open reference systems, respectively. The 
fact that real salt concentrations are generally lower 
in the suspension than in the reservoir [65] indicates 
a further unphysical property of the open reference 
system. 

Accuracy of the linearization approximation can 
be probed by self-consistency checks. First, the av- 
erage linearized potential in the suspension can be 
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FIG. 3: Linear-response theory predictions for osmotic 
pressure (a) and phase diagram (b) of aqueous suspen- 
sions of macroions of diameter a = 50 nm and fixed 
effective valence Z = 500 in Donnan equilbrium with a 
1:1 microion reservoir at salt activity z s = 0.1 mM. Per- 
turbation about closed reference system predicts phase 
instability [(a), solid curve] and vapor-liquid coexis- 
tence [(b), solid curve]. Perturbation about open ref- 
erence system (dashed curves) predicts increased stabil- 
ity with increasing reservoir-to-suspension volume ratio 
[V r /V = 0, 0.2, 0.5, 1, 2, 5, oo, bottom to top in (a)] and, 
correspondingly, a shrinking binodal [V r /V = 0, 0.2, 0.5, 
largest to smallest area in (b)]. Instability vanishes and 
binodal collapses for V r /V > 1. Circles denote critical 
points and tie lines join coexisting phases. Inset of (a): 
salt concentration of suspension c s with closed (solid) 
and open (dashed) reference systems for infinite reser- 
voir (Vr/V — OO). 



estimated from Eq. (35): 

/3e|*| = 47r/3e^ / drr 2 |*(r)| 



3 ZXb t) 
(kq) 2 a l — i] 



(50) 




For the electrostatic coupling considered here 
(ZXb/cl — 7), this estimate yields /9e|^| < 0.3 along 



c 1.1 



0.05 0.1 

Colloid Volume Fraction n 

FIG. 4: Microion density inequality [41], n+n_ > n r , 
tested (and confirmed) at same parameters as Fig. 3 and 
salt activities z s — 0.1,0.2,0.4 mM (top to bottom). 



the binodal of Fig. 3(b). Although the potential is 
relatively high at the surface, /3e|^(a)| ~ 3, it decays 
rapidly with distance from the surface. Integration 
of the microion density profiles [66] for these param- 
eters shows that only 10-20% of the counterions are 
in the highly nonlinear region, /?e|\&(r) — > 1, 
for salt activities z s > 0.05 mM. By comparison, 
the artifacts of linearized PB theory demonstrated 
in Refs. [39-42] occur at much stronger couplings 
(ZXb /a > 18 and (3e\^\ > 1), where nonlinear ef- 
fects are more significant. Second, when compared 
with (salt-free) primitive model simulations [18], 
linear-response theory [64] accurately predicts os- 
motic pressures for the same strength of electrostatic 
coupling, ZX B /a = 7.1 (Z = 40,A B /a = 0.1779). 
Third, the inclusion of first-order nonlinear terms 
in the effective interactions has been shown [35] for 
these parameters to only quantitatively alter the 
phase diagram, although higher-order nonlinear cor- 
rections may not necessarily be small. Fourth, the 
difference in chemical potentials between microion 
species [from Eq. (49)] is found to be negligible. Fi- 
nally, a rigorous inequality involving the average mi- 
croion densities in the suspension [41], n+n_ > n^, 
which follows from Eq. (24), is safely satisfied, as 
demonstrated in Fig. 4. Nevertheless, it is conceiv- 
able that charge renormalization could substantially 
modify the predicted phase behavior, even for such 
a relatively weak coupling. 

Returning to the choice of reference system, the 
common pcrturbative origin of linearized PB and 
response theories suggests a natural selection crite- 
rion. To minimize the magnitude of the perturbation 
term in the microion grand potential [Eq. (3)], and 
thereby optimize the accuracy of linearization, the 
reference system alone should reasonably describe 
the real suspension. On this basis, global electroneu- 
trality clearly favors the closed over the open ref- 
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erence system for a bulk suspension. To appreci- 
ate why, consider that even monovalent macroions 
(Z = 1) would generate a macroscopic charge, which 
must be neutralized by a compensating counterion 
charge, lest the energy density diverge. Only a ref- 
erence system whose counterions are confined to the 
same volume as the macroions, with a density hence 
slaved to the macroions, respects this constraint. 
Electroneutrality thus dictates the reference system 
of a bulk suspension and decouples phase behavior 
from boundary conditions. 

A caveat to the above conclusions may be re- 
quired in the case of Donnan equilibrium of lower- 
dimensional colloidal assemblies, e.g., monolayers 
(quasi-2D) or clusters (quasi-OD), whose ions are all 
near the surface and in direct contact with a reser- 
voir. In such systems, electroneutrality does not 
necessarily constrain the majority of counterions to 
the volume occupied by the macroions. The distri- 
bution of counterions between the macroion region 
and the reservoir depends instead on the strength of 
macroion-counterion attraction and may determine 
whether a suspension is best modeled by an open or 
closed reference system. 

For a colloidal monolayer (e.g., at an air- water 
interface or a glass surface), the range of coun- 
terion localization, compared with layer thickness, 
can be roughly quantified by estimating the Gouy- 
Chapman length I, i.e., the typical distance to which 
thermal energy can separate a counterion from the 
monolayer surface by working against electrostatic 
attraction of the macroions. Approximating a mono- 
layer as a uniform surface charge density, and ne- 
glecting electrolyte screening, yields a lower limit of 
I ~ {i{olZ\bp)~ x , where a is the fraction of escaped 
counterions and p is the areal density. In the dilute 
limit (pa 2 <C 1), the Gouy-Chapman length can far 
exceed the monolayer thickness (I 3> a), implying 
weak localization of counterions by the macroions. 
In contrast, a suspension that is macroscopic in all 
three dimensions, always confines the vast majority 
of its counterions in bulk, considering that I is a mi- 
croscopic length. 

For a small cluster of macroions, a natural mea- 
sure of counterion confinement is the escape energy, 
i.e., the energy required to remove a counterion from 
the cluster to the surrounding reservoir. An iso- 
lated, spherical cluster of N m macroions with total 
charge aN m Ze and volume fraction r\ binds a coun- 
terion near its surface with an escape energy E esc 

that has an upper limit (3E esc ~ a(Z\s /a)Nm 3 r] 1 / s , 
again neglecting electrolyte screening. While coun- 
terions are tightly bound by large, densely charged 
clusters, they can escape from sufficiently small, di- 
lutcly charged clusters via thermal evaporation into 



the reservoir. 

These considerations suggest that the effective in- 
teractions governing lower-dimensional colloidal as- 
semblies may arise from a reference system depen- 
dent on particle number and charge density. A 
closed reference system - appropriate at high N m , 
77, and Z - clearly entails density-dependent effective 
interactions and many-macroion cohesion in deion- 
ized suspensions. In contrast, an open reference sys- 
tem - appropriate at lower N m , ry, and Z , where the 
counterions are not enslaved by electroneutrality to 
the macroions - may yield density-independent ef- 
fective interactions and pair repulsion of macroions. 

Some support for this view comes from observa- 
tions of quasi-2D monolayers in deionized aqueous 
suspensions [68-70]. In these experiments, highly 
charged PS particles were confined to a plane by 
light pressure and corralled via optical traps to ef- 
fective surface charge densities of 0(1O -5 ) C/m 2 or 
0(1O -4 ) e/nm 2 . Effective pair potentials, obtained 
by inverting measured radial distribution functions, 
exhibited a density-independent Yukawa form at 
lower densities, crossing over to density-dependent 
non- Yukawa behavior at higher densities. These 
observations are consistent with salt-dominated lin- 
ear screening in dilute monolayers, whose microion 
population is fixed by the reservoir, and nonlin- 
ear screening above a threshold macroion density. 
Colloidal monolayers confined between glass plates 
also reportedly can exhibit density-dependent effec- 
tive pair interactions [71], although interpretation 
of such experiments is complicated by uncontrolled 
microion exchange with the glass surfaces. 

Some evidence for size dependence in the sta- 
bility of colloidal clusters comes from observations 
of metastable crystallites and macroion gathering 
in low-ionic-strength suspensions compressed near 
glass plates by external electric fields [7, 8]. Slow 
counterion evaporation from the macroions, and the 
associated crossover from cohesive many-body to 
repulsive pair effective interactions, might qualita- 
tively explain the slow break-up of small macroion 
aggregates adrift in deionized water. For example, 
the fee crystallites reported in Ref. [7] - character- 
ized by a = 326 nm, Z = 7300, \ B = 0.72 nm, 
c s < 5 jtiM, and nearest-neighbor separations d = 

2 /3 

1.8 /im - would correspond to [3E CSC < 0.56aN„{ . 
A crystallite of size N m — O(10 2 ) could conceivably 
break up by gradually losing a substantial fraction of 
its counterions to the surrounding reservoir. Such a 
destabilization mechanism also could account for the 
observed increase in stability with increasing cluster 
size, without invoking any long-range pairwise at- 
tractive interaction. Although a crossover between 
closed and open reference systems appears consis- 
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tent with these experimental observations, a quan- 
titative description requires further work. 

V. CONCLUSIONS 

Summarizing, this work analyzes the influence of 
the reference system on the effective electrostatic 
interactions and phase behavior predicted by lin- 
earized, coarse-grained theories of charged colloids. 
For a suspension in Donnan equilibrium with a mi- 
croion reservoir, two reference systems are conceiv- 
able: one closed and the other open with respect 
to microion exchange with the reservoir. The con- 
straint of global electroncutrality is shown to pro- 
vide an objective physical basis for selecting between 
these two reference systems. From this criterion, and 
the unification of Poisson-Boltzmann and response 
theories, it is concluded that bulk suspensions are 
properly modeled within PB theory by expanding 
about the average potential and average microion 
densities within the suspension (not the reservoir), 
and within response theory by perturbing about a 
closed reference system, whose countcrions share the 
same volume as the macroions. 

Dependence of predictions on the reference system 
need not necessarily constitute a failing of linearized 
theories; it merely underscores the importance of 
correctly choosing the reference system. The op- 
timal reference system identified in Ref. [41] for PB 
theory is the same one that respects electroneutral- 
ity in linear-response theory. The present conclu- 
sions thus broaden and clarify those of previous stud- 
ies based on PB cell models [39-42] . When linearized 
about a physically consistent reference system, both 



PB and response theories predict a spinodal phase 
instability of highly dcionized suspensions, assum- 
ing a state-independent effective macroion charge. 
Whether the predicted instability has any relevance 
to experimentally observed anomalies, however, can 
be decided only by resolving the important issues of 
nonlinear screening and charge renormalization. 

Lower-dimensional systems - from monolayers to 
clusters - may be best described, depending on the 
strength of macroion-counterion correlations, by a 
closed reference system at high charge 1 densities, but 
an open reference system at lower charge densities. 
The implied transition in effective interactions and 
thermodynamic properties may parallel an observed 
crossover, as a function of density, in the structure of 
colloidal monolayers and might qualitatively explain 
reports of metastable crystallites. Further compar- 
isons with experiment are required to more fully 
evaluate these speculations. Future work should 
seek a more unified, cross-dimensional description 
of charged colloids. Generalization of the coarse- 
grained approach to richer models with multiple mi- 
croion species - some exchanging with a reservoir 
and some trapped - may have applications to col- 
loids or nanoparticles confined in pores and to pro- 
teins inside cells, where ion channels regulate ex- 
change of ions across cell membranes. 
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